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ABSTRACT 

We report the detection of infrared emission from the jet of the nearby FR I radio galaxy 3C 31. 

The jet was detected with the IRAC instrument on Spitzer at 4.5 /rm, 5.8 /im, and 8.0 /im out to 
30" (13 kpc) from the nucleus. We measure radio, infrared, optical, and X-ray fluxes in three regions 
along the jet determined by the infrared and X-ray morphology. Radio through X-ray spectra in 
these regions demonstrate that the emission can be interpreted as synchrotron emission from a broken 
power-law distribution of electron energies. We find significant differences in the high energy spectra 
with increasing distance from the nucleus. Specifically, the high energy slope increases from 0.86 to 
1.72 from 1 kpc to 12 kpc along the jet, and the spectral break likewise increases in frequency along 
the jet from 10-100’s of GHz to ~20 THz. Thus the ratio of IR to X-ray flux in the jet increases by at 
least an order of magnitude with increasing distance from the nucleus. We argue that these changes 
cannot simply be the result of spectral aging and that there is ongoing particle acceleration through 
this region of the jet. The effects of mass loading, turbulence, and jet deceleration, however these 
processes modify the jet flow in detail, must be causing a change in the electron energy distribution 
and the efficiency of particle acceleration. 

Subject headings: galaxies: individual(NGC 383)— galaxies: active — galaxies: jet — infrared: galax¬ 
ies 


1. INTRODUCTION 

Extended radio emission in galaxies is the observable 
manifestation of powerful, collimated outflows from su- 
permassive black holes. Only ~10% of active galaxies are 
radio loud and contain well-formed jets, but in all but the 
most luminous systems, the mechanical power of the jets 
gene rally dwarfs any rad iative losses from an accretion 
disk (|Birzan et al. |[2004ll . It is now well established that 
jets play an important role in feedback between the cen¬ 
tral supermassive black hole (SMBH) and its larger scale 
environment. However, in spite of several decades of re¬ 
search, both observational and theoretical, basic issues 
related to jet physics, including particle content, veloc¬ 
ity structure, the role of magnetic fields, and a host of 
other problems, remain unresolved. Extragalactic jets 
in low-power Fanaroff-Riley type I (FR I) sources are 
generally only observable at frequencies above the radio 
band via the synchrotron radiation from ultra-relativistic 
particles accelerated at shocks or other disturbances in 
the flow, although at X-ray energies, inverse Compton 
emission can sometimes be important. The significant 
numbers of detections of FR I jets in the IR, optical, and 
X-ray regimes has given us new insig hts into jet dy nam- 
ics and particle acceleration (see e.g. lWorrallir2009l for a 
review). 

The radio, optical, and X-ray emission from FR I jets 
can normally be modeled as synchrotron emission from 
a single electron energy distribution described by a bro¬ 
ken power law. Electrons emitting synchrotron radia¬ 
tion have lifetimes that vary inversely with their energy, 
so multi-wavelength observations are expected to show 
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that higher-energy electrons are concentrated closer to 
the sites of particle acceleration than their lower-energy 
counterparts. In typical equipartition magnetic fields of 
FR I radio galaxies, the lifetimes of IR/optical and X- 
ray emitting partic les are ^10^ and 10^ yrs, respectively 
('e.g. lWorrallll2009ll . Light travel times for kiloparsec jets 
demonstrate that the ultra-relativistic particles must be 
accelerated in situ. Proposed acceleration mechanisms 
include magnetohydrodynami c turbulence a nd shock ac¬ 
celeration, among others le.g. lWorrallll2009l l. 

3C 31 is a nearby twin-jet radio source hosted within 
the early-type galaxy NGC 383, the dominant member 
of a small group of galaxies that also includes the nearby 
companion NGC 382. At a redshift of 0.0168 (11^=72.9 
Mpc) (|de Vaucouleurs et al. IHOQUI^ . 3C 31 was one of 
the ea rliest FR I galaxies to be studied in depth in the 
radio (jBnrchl 119771 iFomalont et al. 1 11980( 1 and remains 
a canonical example of this class. The twin jets are 
aligned roughly along a north-south axis. The north¬ 
ern jet is signific antly brighter at G Hz frequencies in the 
central 10 kpc (|Lara et al. I Il997ll . suggesting, that in 
this regime at least, the flow is mo derately relativisti c 
and this jet is oriented towards us (iLaing et al. II200§1 . 
The northern jet is straight and narrow within 45" (19 
kpc) of the nucleus, but beyond this distance, flares 
dramatically suggesting that it is losing kinetic energy 
(|Burchlll977f) . On larger scales, both the forward and re- 
yerse jets extend at least 20' (500 kpc) from the nucleus 
(|Laing et al. 1 1200§1 . The radio morphology of the jet 
becomes wavy and diffuse at these distances. Subsonic 
motions of the external medium are probably important 

® We adopt a cosmology with Hq = 70 km Mpc~^, Qm = 
0.27, and Ua = 0.73, resulting in a linear projected scale of 0.34 
kpc per arcsec at the source. W e deproject distances as suming the 
jet is 52° from the line of sight IILaing fc Bridld[2002bl'l . 
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in de termining the large scale radio morphology of the jet 
le.g. lLoken et al. II1995D . The inner 8.5 k pc of the north¬ 
ern jet has been detected in the optical (|Croston et al. I 
I2008L[pi tcher et al. 1119801 ). and the inner 2.7 kpc in the 
X-ray ( Hardcastle et al. 1 1200211 . It has also previously 
been detected a t lower spatial resolution in the infrared 
(|Tanslevlll997tl . iHardcastle et al. I (|2002li modeled the jet 
emission as synchrotron radiation from a population of 
ultra-relativistic electrons. 

The relatively regular two-sided nature of 3C 31 in 
radio has permitted detailed kinematic modeling of the 
jets. In particular, Laing and Bridle (2002a, 2002b, 2004) 
created a 2-dimensional model of the jet flow to deter¬ 
mine the velocity profile as a function of distance down 
the jet axis. Laing and Bridle (2002a, 2004) extended 
their work to include dynamical modeling by incorporat¬ 
ing the pressure of the exterior gas as a function of dis¬ 
tance based on X-ray measurements. They divided the 
inner jet into three regions: an inner region (extending 
to 2'.'5 (1.1 kpc)) wherein the jet expands as a cone with 
a half-opening angle of 8?5, a flaring region (2'.'5 — 8'.'3; 
1.1 - 3.6 kpc) with gradually larger opening angle, and 
an outer region (8'.'3 — 28'.'3; 3.6 - 12.2 kpc) wherein the 
jet once more expands conically, but with a half-opening 
angle of 16? 75. They found that these regions have dis¬ 
tinct kinematic properties and that the jet deceleration 
is consistent with the entrainment of gas from the X-ray 
emitting medium. 

In this paper, we present a detection of infrared emis¬ 
sion from the jet of 3C 31 as imaged by the Spitzer Space 
Telescope Infrared Array Camera (IRAC). Spitzer has 


detected infrared lets in a few nearby galaxies, includ- 

ing M87 (IShi et al. 1 120071. iForman et al. 

20071! and Cen 

A (NGC 5128! (IHardcastle et al. 1120061 

Brookes et a,l. 1 


HOOI). We discuss 3C 31’s inner jet morphology in the 
infrared. X-ray, and radio. We measure jet IR flux den¬ 
sities in the IRAC and MIPS bands in the fla ring and 
outer regions defined bv iLaing fc Bridld (l2002al) . We use 
the resulting energy distribution for each region to con¬ 
strain particle acceleration in the jet. Images have North 
to the top and East to the left, angles are given counter¬ 
clockwise from North, and spectral indices a are defined 
such that flux density at frequency v oc v~°^. We use 
th e 52° angle to the line o f sight for the jet determined 
bv iLaing fc Bridl^ (l2002bll (hereafter L02) to give linear 
distances as deprojected distances. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Spitzer Data 

IRAC (|Fazio et al. Il2004) observations of 3C 31 were 
obtained on 2006 February 15 and 22 (Spitzer PID 3418) 
at 3.6 /rm, 4.5 /xm, 5.8 /xm, and 8.0 ^m. Each obser¬ 
vation consisted of 24 frames of 30 s integrations cov¬ 
ering the galaxy and the nearby field obtained using a 
medium cycling dither. For our analysis, we retrieved 
the Basic Calibrated Data (BCD) version S14.0 pipeline 
products from the Spitzer archive. For the 8.0 /xm obser¬ 
vation, light scattered from the galaxy nucleus along the 
detector-array rows was fit and subt racted from the BCD 
frames using custom software as in lAshbv et al. I (I2009D . 
We then coad ded the 30 s BCD fram es using version 4.1.2 
of IRACproc (|Schuster et al. 11200611 into mosaics having 
1"2 pixels, i.e., the native IRAC pixel size. A 4.'8 x 4.'8 


region centered on the galaxy was selected for analysis 
in each IRAC mosaic. 

To detect and determine the jet’s infrared morphol¬ 
ogy, we needed to subtract the bright infrared stellar 
emission. We modeled the stellar emission using the 
3.6 / xm emission, as has been done for o ther galaxies 
fe.g. iPahre et aZ. 112004 iTang et a/.'1l2009D . After sub¬ 
tracting a constant background determined in the frames 
far from NGC 38 3 and NGC 382, we used GALFIT 
(|Peng et al Il2002h to fit the ellipticity and position an¬ 
gle of the galaxy. The 3.6 ^m image was scaled to cor¬ 
rect for differing instrumental parameters between the 
different IRAC bands, specifically different aperture cor¬ 
rection factors for an infinite aperture and the different 
photometric zeropointfl We then performed photome¬ 
try using elliptical annuli on the 4.5 ^m, 5.8 ^m, and 
8.0 /xm background-subtracted images and the associ¬ 
ated scaled 3.6 /xm images. We measured the flux in 
each annulus and used these to calculate the colors be¬ 
tween 3.6 /xm and each of the longer IRAC bands in each 
annulus, resulting in color profiles for NGC 383. We fit 
these profiles in the annuli from 6'.'5 to 11875 centered on 
the nucleus. When extrapolated to the galaxy core, the 
fits were used to color-correct the scaled 3.6 /xm images 
to be the stellar models at 4.5 /xm, 5.8 /xm, and 8.0 /xm. 
This method does not take into account differences in 
colors due to central features, so we expect that some 
residual emission may be pre sent from eith e r the AGN 
or the dust disk described bv iMartel et a.l. I ((20001! . We 
detected nonstellar jet emission in the 4.5 /xm, 5.8 /xm, 
and 8.0 /xm bands, after the stellar model was subtracted. 
The residual images at 5.8 /xm and 8.0 /xm are shown in 
Figure [TJ The 8.0 /xm image is also shown in Figure 
[2] with contours of radio emission in order to localize 
the infrared emission within the larger radio jet. Fig¬ 
ure [3] shows the 4.5 /xm residual image along with a 
3.6 /xm image, where an azimuthally smoothed version of 
that image has been subtracted in order to detect the jet. 
While the jet is less clearly present at 4.5 /xm than in the 
two longer IRAC bands, these residual images in the two 
shorter IRAC bands contain emission at the location of 
the jet. The technique used to subtract the stellar emis¬ 
sion of NGC 383 results in different residuals for NGC 
382, but does not impact our study of the jet. We also 
modeled the stellar emission using 2MASS K-band emis¬ 
sion. Finally, we used the IRAF ELLIPSE package to 
model the galaxy emission (Bliss et al., in preparation). 
The resulting nonstellar images show similar morphology 
to those in Figures [T]131 

We also obtaine d a mosaic at 24 /xm taken with the 
MIPS instrument (|Rieke et al. 11200^ on 2004 December 
26 (Spitzer PID 82). We used GALFIT to fit the emission 
with a galaxy model along with two point sources for the 
nuclear emission of 3G 31 and for the companion galaxy 
NGC 382. 


2.2. Chandra, GALEX, Hubble, 2MASS, and VLA Data 

3C 31 was observed for 44.41 ks on 2000 November 6 
(ObsID 2147) with th e Chandra Advanced CC D Imaging 


Spectrometer (ACIS: IWeisskoof et al 

120001!. This ob- 

servation was previously analyzed by 

Hardcastle et al. 1 


^ IRAC Instrument Handbook Tables 4.1 and 4.7: 
http: //ssc.spitzer.caltech.edu/irac/iracinstrumenthandbook/ 
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Fig. 1.— Spitzer IRAC images of the nonstellar emission at 5.8 /^.m (left) and 8.0 fim (right) in a lf5 x 1' field around the nucleus. The 
stellar emission was modeled using the 3.6 //m scaled to compensate for zero point offsets and aperture corrections. The 3.6 /im emission was 
color-corrected using a color profile determined for each of the longer IRAC bands. Jet emission is seen in the region to the north-northwest 
of the nucleus, both of which are indicated by the red arrows. The technique used to subtract the stellar emission of NGC 383 results in 
different residuals of NGC 382 but does not impact our study of NGC383’s jet. The scale bar is 15^^ (6.5 deprojected kpc). 



Fig. 2.— Spitzer 8.0 /im image of the jet of 3C 31 in 1^2 pixels with 8.4 GHz radio contours overlaid, showing the relative location of 
the IR emission in the larger radio jet. The image was smoothed with a Gaussian 0('6 wide after the stellar emission was subtracted. The 
8.4 GHz contours are at 0.1, 0.2, 0.4, 0.6, 1.2, 2, 20 mjy beam”^ 


IMl). GALEX observed 3C 31 on 2006 November 9 
(AIS 42) for 109 se conds in both the NUV and FUV 
filter ([Martin et al. I [20051) . We obtained the intensity 
maps in these filters from the Multimission Archive 
at STScI (MAST). We retrieved mosaics of 3C 31 
taken by the Wide Fiel d Planetary Camera 2 (WFPC2) 
(jHoltzman et al. 1119951 ) on the Hubble Space Telescope 
(HST) from the Hubble Legacy archives. 3C 31 was ob¬ 
served through the F555W and F814W filters on the 
PCI chip for 460 seconds each on 1998 September 23 
and through the FR680N filter on the WF2 chip for 600 


(a' 69 X 0"67 beam). 

seconds on 1995 September 01. 2MASS All-Sky Survey 
Atlas images through all three filters were obtained from 
the NASA/IPAC Infrared Science Archive (IRSA)@. The 
VLA data we use here are observations for 6.5 hrs at 
8.4 GHz on 1994 June 14 in the B configuration and for 
25.5 minutes at 1.4 GHz on 1999 July 18 in the A con¬ 
figuration. We used the Astronomical Image Processing 
System (AIPS; version 31DEC09) to generate maps with 
beam sizes of 0"69 x 0"67 at 8.4 GHz and 1'.'76 x 1'.'40 at 

® http://irsa.ipac.caltech.edu/applications/2MASS/IM/ 
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1.4 GHz. 

3. RESULTS 
3.1. Jet Morphology 

The infrared jet emission detected at 5.8 /rm and 
8.0 /im is shown in Figure [TJ The jet and nucleus are 
identified with red arrows. The northern jet is detected 
from 8" (3.5 kpc) to 24" (10.3 kpc), with a position an¬ 
gle of —25°. The X-ray emission from the northern radio 
j et ha s previously been described by iHardcastle et al. I 
(|2002tl . It extends approximately 8" (3.5 kpc) from the 
nucleus. We have a 3.5cr detect ion of X-ray e mission from 
8"5 (3.5 kpc) to 15" (6.5 kpc). iBurchl (I1977D mapped 3C 
31 at 5 radio fr equencies between 0 .4 GHz and 15 GHz. 
More recently, iLainsr et al. I (j2008ll obtained high reso¬ 
lution VLA images from 1365-8440 MHz. An infrared 
image of the jet with 8.4 GHz contours of the inner ra¬ 
dio emission is shown in Figure [5J The radio emission 
shows a pair of jets, whose width increases with distance 
from the nucleus. The northern jet is brighter than the 
southern jet by about an order of magnitude between 
2'.'5 — 6" (1.1-2.6 kpc) and remains straight for ~45" be¬ 
fore gently bending to the north, while the southern jet 
extends 120" before exhibiting a strong bend. The jet 
is not detected in th e Hubble or GALEX observations. 
iGroston et al. I (I2003D detected jet emission using the 4 
m Mayall Telescope in the B and R bands. Images of the 
jet emission in the radio, IR, and X-ray, as well as the 
corresponding region in the optical (with HST) and UV 
(with GALEX), are shown in Figure 01 

To quantify the significance of the detection of the 
northern jet in the infrared, we measured the azimuthal 
surface brightness profile in an elliptical annular sector 
with an ellipticity of 0.87 and with inner and outer semi¬ 
major axes of 12" (5.2 kpc) and 20" (8.6 kpc), where the 
peak of the infrared jet emission is located based on the 
residual images shown in Figure [TJ The total 8.0 /rm sur¬ 
face brightness in this annulus is shown in Figure [5] The 
statistically significant surface brightness enhancement 
between 330° and 350° is due to the northern jet emis¬ 
sion. 

In contrast to the northern jet, the counterjet is only 
marginally detected in the two longest wavelength IRAC 
bands. Faint IR emission is noted in the Spitzer image 
shown in Figure [21 but the statistical significance of this 
emission is low (less than 3cr). The point estimate of the 
flux density of the counterjet is sensitive to the details 
of the model used to subtract the stellar emission, but 
the statistical significance of the detection is insensitive 
to the stellar model. We measured the flux in the IRAG 
bands in a 10" x 5" (4.3 kpcx2.2 kpc) aperture begin¬ 
ning 15" from the nucleus centered in the radio contours. 
We find flux densities with Icr uncertainties of 35 ± 15 
/rjy at 5.8 /rm and 16 ± 16 /rJy at 8.0 ^m. In compar¬ 
ison, the forward jet flux density at the same distance 
from the nucleus is significantly greater at 85 ± 16 /rJy 
at 5.8 /im and 142 ± 10 /iJy at 8.0 /im. 

Figure [6] shows the longitudinal profile of the northern 
jet emission at 8.4 GHz, 8.0 /im, and in the soft X-ray 
band (0.5-2.0 keV). In measuring the X-ray profile, we 
removed the emission from the ambient hot gas by sub¬ 
tracting a profile in similar rectangles adjacent to the jet. 
We do not try to deconvolve the Spitzer emission within 


7" of the nucleu s, which cont a ins th e dusty disk seen 
in the optical by iMartel et al. I (|2000[1 . as well as stellar 
emission and an unresolved point source. The profiles 
in Figure [6] were measured in rectangles 4" wide and ei¬ 
ther 1'.'2 (IR and radio) or 2" (X-ray) long along the 
jet axis. The different bands have different beam sizes, 
which results in an oversampled profile for the Spitzer 
8.0 fj,m band. We binned the data in fewer regions in the 
X-ray band to emphasize the appearance of this emis¬ 
sion. The 8.0 /im profile shows that the infrared jet is 
brightest bet ween 12" (5.2 kpc) an d 20" (8.6 kpc) from 
the nucleus. IGroston et al. I (|2003l l measured the longi¬ 
tudinal profile of the jet in the B band (their Figure 5). 
They detected optical emission from the jet within ~ 25" 
of the nucleus with an increase in brightness coincident 
with the brightest infrared jet emission that we observe. 
The radio profile shows a slow decrease in brightness with 
distance along the northern radio jet, with a bright re¬ 
gion where the X-ray emission is strongest closest to the 
nucleus. 

3.2. Jet Spectrum and Electron Energy Distribution 

We measured the radio, infrared, and X-ray flux den¬ 
sities from the jet in three different regions (labeled Rl, 
R2, and R3; see Figure [7]) along the jet axis and set upper 
limits on the UV and optical fluxes in these regions. Our 
choice of regions was guided by the morphology of the 
jet in the IR and by its relationship to the X-r ay emis¬ 
sion seen with Ghandra (|Hardcastle et ui. ]l2002fl . In the 
dynamical study of the 3G 31 jet by L02, the jet was 
divided into three regions based on its radio morphol¬ 
ogy and modeled jet dynamics. These three regions were 
termed the inner, the flaring, and the outer regions, and 
differ somewhat from our regions. Region Rl contains 
the peak X-ray emission and extends from 275 (1.1 kpc) 
to 873 (3.6 kpc) from the nucleus with a width of 4" 
(1.7 kpc). This roughly corresponds to the flaring region 
in the analysis of L02. We divided the outer region as 
defined by L02 into two distinct regions, which we call 
R2 and R3. The R2 region extends from 873 (3.6 kpc) 
to 1873 (7.9 kpc) from the nucleus, and the R3 region 
extends from 1873 (7.9 kpc) to 2873 (12.2 kpc). Both 
regions are 5" (2.2 kpc) in width. The inner region of 
L02, where the jet has roughly constant opening angle, 
is within 275 of the nucleus and is not resolved from the 
host galaxy in the Spitzer data, so we did not use it in 
our analysis. 

In addition to the radio, IR, and X-ray detections of the 
jet, UV, optical, and additional IR data points are taken 
from archival GALEX, HST, and 2MASS and MIPS ob¬ 
servations, respectively. The fluxes with uncertainties 
and upper limits at all frequencies for the three regions 
are given in Table [T] We measured consistent fluxes 
in the IRAG bands from the nonstellar images result¬ 
ing from the subtraction of the IRAF ELLIPSE galaxy 
model and of the 3.6 /rm stellar emission model. 

The broadband spectra from Rl, R2, and R3 are shown 
in Figures [8] through [101 To better constrain the elec¬ 
tron energy distributions (EEDs) in the three regions and 
quantify differences between them, we fit the radio to X- 
ray spectra with power-law models. In all the model fits, 
we used only the actual detections with the constraint 
that the ht model could not exceed any of the upper 
limits. The limited number of points in the spectra and 
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Fig. 3.— Spitzer IRAC images at 3.6 fim (left) and 4.5 fim (right) in a square arcminute field around the nucleus. The 3.6 fim image 
has had an azimuthally smoothed version of the image subtracted, in order to emphasize the jet emission. The 4.5 pm image is of the 
nonstellar emission where the stellar emission was modeled based on the scaled 3.6 pm emission in the same manner as for Figure[T] 



Fig. 4.— Jet emission at 1.4 GHz (a), 8.4 GHz (b), nonstellar 8.0 pm (c), and broadband X-ray energies (0.3-8.0 keV) (f). Optical 
(FR680N, WFPG2 on HST; shown in log scale) and NUV (GALEX) images of the same region are shown in (d) and (e) respectively. The 
15” scale corresponds to 6.5 kpc. 


the relative simplicity of the models does not warrant a 
more sophisticated statistical treatment of the upper lim¬ 
its. Only some of the upper limits put strong constraints 
on the model fits. We first fit single power-law models to 
the three spectra and found that this simple model could 
be rejected at >99.9% significance in all cases. 

We therefore fit the spectra with models of synchrotron 
emission from broken power-law electron energy distri¬ 
butions which contain four free parameters: two spectral 
indices, a break frequency, and a normalization. An ac¬ 
ceptable fit is found for the spectra from all three regions 
with this model. The model is overdetermined in the R1 
region since we have only three data points and a model 
with four free parameters. In this region, we consider 
two models with the break frequency frozen. The optical 
upper limit from HST forces the break frequency to lie 
in the optical or at lower frequencies, so we consider two 


models that bracket the extreme limits - one model (FI) 
with a break frequency just above the 8.4 GHz radio data 
point and a second model (F2) with a break frequency 
in the mid-IR, the highest frequency that fits through 
the radio and X-ray detections and lies below the opti¬ 
cal upper limit. The best fit parameters for each of our 
three regions are summarized in Table [2l and the best 
fit emission models have been overplotted on Figures [5] 
to [TUI 

There are a few interesting points to note about the 
morphology of the jet in the radio, IR, and X-ray bands. 
X-ray emission is detected only in the R1 and R2 regions, 
the two regions closest to the nucleus. The presence of 
the X-ray emission suggests that there is on-going par¬ 
ticle acceleration in these regions, but that acceleration 
to the highest energies appears to be less efficient in the 
outer parts of the jet. Guided by our experience with 
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TABLE 1 
Jet Fluxes 


Wavelength 

R1 

(^Jy) 

Flux Density in Region 
R2 

(nJy) 

R3 

(iiJy) 

1.4 GHzi 

7.87 X lO'^ 

8.98 X lO'^ 

6.73 X 10-* 

8.4 GHzi 

3.01 X lO'^ 

3.85 X 10^ 

2.63 X 10-* 

24 fim 

< 960 

< 709 

< 651 

8.0 fim 

< 1200 

152 ± 25 

137 ±9 

5.8 fim 

< 800 

99 ± 20. 

78± 19 

4.5 fim 

< 700 

< 71 

50. ±7 

3.6 fim 

< 1200 

< 94 

< 48 

2.2 fim 

< 2400 

< 570 

< 260 

1.7 fim 

< 2500 

< 510 

< 270 

1.2 fim 

< 1800 

< 570 

< 380 

820.3 nm 

< 6.9 

< 10. 

< 10. 

667.3 nm 

< 17 

< 25 

< 25 

520.2 nm 

< 4.8 

< 6.9 

< 6.9 

226.7 nm 

< 12 

< 10. 

< 7.5 

151.6 nm 

< 19 

< 10. 

< 8.4 

0.3-8.0 keV 

(3.3 ±0.3) X 10 

-3 (3.1 ± 1.0) X lO-'^ 

< 1.3 X lO-'* 


Note. — The R1 region extends 2^5 — 8^3 and is 4'' wide. The R2 region 
extends — 18^3 and is b" wide. The R3 region extends 18^3 — 28^^3 and 
is 'o" wide. The upper limits are 3cr upper limits. 

^ The formal statistical uncertainties are smaller than the significant digits. 
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Fig. 5.— Azimuthal profile of the total 8.0 //m emission for an 
elliptical annulus with inner and outer radii of 12" (5.2 kpc) and 
20" (8.6 kpc) with an ellipticity of 0.84. Angles are given counter¬ 
clockwise from north. The jet is located at 330° — 350°. 


Cen A (|Kraft et al. 112001 IHardcastle et al. II2003I) . there 
is almost certainly unresolved substructure in these two 
regions. This clearly demonstrates that the electron en¬ 
ergy distributions, and probably the nature of the parti¬ 
cle acceleration in general, are distinctly different in each 
of these regions. There are significant differences in the 
three spectra that clearly define three distinct regions 
of particle acceleration. The low-frequency (i.e. radio) 
spectral indices for the three regions are between 0.47 
and 0.54, typical values for radio plasma. The break 
frequency in region R1 is poorly defined, but must lie 
somewhere between 10-800 GHz. The spectral index af¬ 
ter the break is relatively flat, ~0.86-0.97. Regions R2 
and R3, in contrast, have fairly well defined break fre¬ 
quencies around 10^^ Hz, in the far-IR spectral regime, 
with a significant spectral steepening beyond the break. 
In particular, the spectral index after the break is ~1.28 
and ^1.72 in regions R2 and R3, respectively. 

Using the best-fit broken power-law models for the 
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Fig. 6.— Longitudinal profile of the northern jet emission at 
8.4 GHz, 8.0 fim and 0.5-2.0 keV. Flux densities were measured in 
rectangles 4" wide and 1'.'2 long along the jet axis for the infrared 
and radio. The 8.0 fim emission was taken from the nonstellar 
image (Fig. X-ray flux densities were measured in rectangles 
4" wide and 2" long. We removed the X-ray emission from hot 
ambient gas, determined in similar rectangles adjacent to the jet. 
The regions that will be used to extract spectra are marked at the 
top of the plot. 

spectra, we determined the equipartition magnetic fields 
and the break energies of the electron distributions. Pre¬ 
dictions for X-ray IC/CMB and SSC e mission have been 
computed using the synch program (iHardcastle et al] 
[HUllOOl. We assume that 7mm=10 and jmax=W^, 
K=0 (the ratio of non-radiating to radiating particles), 
that the radio plasma has unity filling factor, and that 
the emitting volumes can be modeled as cylinders in the 
plane of the sky. The corrections for bulk relativistic 
motion were small, even in R1 where the jet velocity 
is greatest, due to the relatively large angle to the line 
of sight. In the two regions where X-ray emission is 
detected, the observed X-ray flux density lies orders of 
magnitude above the IC/CMB and SSC predictions, thus 
strengthening the conclusion of lHardcastle et alTI (|2002f ) 
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Fig. 7.— The apertures used for measuring the spectra overlaid on the 1.4 GHz (a), nonstellar 8.0 /im (b), and broadband X-ray (0.3- 
8.0 keV) (c) emission. Note that there is radio emission in all three apertures and that the IR jet emission is concentrated in the two 
regions farthest from the nucleus, while the X-ray emission is predominantly in the region closest to the nucleus. The 15'' scale corresponds 
to 6.5 kpc. 


TABLE 2 

Broken Power Law Fit Parameters 


Region 

Normalization 

(Jy) 

a 

Low Frequency 

a 

High Frequency 

Break Frequency 
(GHz) 

Rl-Fli 

6300 

0.54 

0.86 

8.42 

RI-F 2 I 

6300 

0.54 

0.97 

780^ 

R2 

19001 

0.471 

1.28 ±0.05 

5250 ± 1330 

R3 

42001 

0.521 

1.72 ±0.12 

18600 ± 1500 


Note. — Fit parameters obtained by minimization using detected fluxes (see Table 
rriand requiring the upper limits to be respected. Since the R1 model is overdetermined, 
we fix the break frequency at the two extremes permitted by the upper limits and de¬ 
termine the remaining parameters. The uncertainties quoted are the 90% uncertainties. 

^ The formal statistical uncertainties of these fits or parameters are negligible compared 
to other sources of uncertainty. 

^ Extremal break frequencies set by the existing radio data and the HST upper limits. 


that this is synchrotron emission from a population of 
ultra-relativistic electrons. We find equipartition mag¬ 
netic fields of 3.0, 2.3, and 2.1 nT in regions Rl, R2, and 
R3, respectively. The equivalent magnetic pressures lie 
somew hat below the pressure of the hot ISM external to 
the jet (|Laing fc Bridlel[2002al) . and implies that there is 
additional pressure support inside the jet, perhaps from 
lower-energy relativistic electrons or relativistic or ther¬ 
mal protons. At these magnetic field strengths, we find 
that the synchrotron lifetime for infrared emitting par¬ 
ticles is ^ 3 X 10^ yr. In order for the observed X-ray 
flux densities to be the result of IC/CMB, the magnetic 
field would need to be two orders of magnitude below its 
equipartition value. 

The observed IR emission from the 3C 31 jet cannot 
be plausibly attributed to any type of inverse-Compton 
or electron scattering process for two reasons. First, if 
the IR emission were not synchrotron radiation, the elec¬ 
tron energy distribution would have to turn down dra¬ 
matically in an equipartition field of 3 nT to explain 
the radio emission but fall below the infrared emission 
(7 ^ 10^“®). The detected X-ray emission in regions Rl 
and R2 could not then be explained by synchrotron emis¬ 
sion from a simple extension of the electron energy dis¬ 
tribution to 7 '^10®. As described above, IC/CMB and 
SSC emission can be ruled out, unless the jet is far from 


equipartition, so the origin of the X-ray emission would 
be unknown. Second, the pressure and total energy of 
thermal electrons from stars in the jet required to create 
the IR emission from electron scattering would be unreal¬ 
istically large. In particular, the energy density of stellar 
photons at the position of the jet would be ^5x10“^"^ 
ergs cm“®, either resulting in an unstable overpressured 
jet or requiring an unrealistically large energy flux due 
to the required number of scatterers. Alternatively, as 
an FR I radio galaxy, 3C 31 is a mis-directed BL La c 
in the Unified Scheme of AGN (lUrrv & PadovanillI995fl . 
There could be a significant component of IR emission 
from the pc-scale jet that is not beamed in our direc¬ 
tion. In this scenario however, a typical BL Lac would 
require an electron density of ^I cm“® in its jet in order 
to account for the observed infrared flux. If 3C 31 were 
a very bright (^2-3 orders of magnitude brighter) mis¬ 
directed BL Lac requiring a lower density in its jet to 
create the observed infrared flux, we would expect other 
morphological features not seen in this system (e.g. the 
IR/optical beamed flux would be visible at large radii 
due to Thomson scattering from the electrons in the hot 
ISM). 


4. DISCUSSION 
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R1 


R3 




log,g(Frequency/Hz) 


log,g(Frequency/Hz) 


Fig. 8 .— Spectrum of the jet in region Rl. The diamonds denote 
the flux densities of the detections, and the arrows the Sa upper 
limits for the non-det ections. The ci r cles a re estimates based on 
the measurements of ICroston et aU\ (I20Q3IV Two representative 
models of the synchrotron output from broken power-law electron 
energy distributions are plotted; one in which the break frequency 
lies just above the 8.4 GHz radio point, and a second in which 
the break frequency has the largest allowable value given the HST 
constraints. The predicted IC/CMB and SSC X-ray flux densities 
assuming equipartition are shown as the dotted and dash-dotted 
lines, respectively. 


R2 



log,g(Frequency/Hz) 

Fig. 9. — Spectrum of the jet in region R2. The diamonds de¬ 
note the flux densities of the detections, and the arrows the 3cr up¬ 
per limits for the non-d etections. The c i rcles are estimates based 
on the measurements of ICroston et aL] H20Q3IV The best fit syn¬ 
chrotron model from a broken power-law electron energy distri¬ 
bution is also plotted. The predicted IC/CMB and SSC X-ray 
flux densities assuming equipartition are shown as the dotted and 
dash-dotted lines, respectively. 


Fig. 10.— Spectrum of the jet in region R3. The diamonds 
denote the flux densities of the detections, and the arrows the 3cr 
upper limits for the non- detections. The circles are estimates based 
on the measurements of ICroston et alTl <120030 . The best fit syn¬ 
chrotron model from a broken power-law electron energy distri¬ 
bution is also plotted. The predicted IC/CMB and SSC X-ray 
flux densities assuming equipartition are shown as the dotted and 
dash-dotted lines, respectively. 


We detected IR emission from two regions of the north¬ 
ern jet in 3C 31 within ~ 15 kpc of the nucleus. The 
third region Rl likely also contains IR emission, how¬ 
ever, we cannot resolve it from the nucleus and dust disk 
also present within this region. From the fitted spectra 
in Rl, we estimate that IR emission is most likely ap¬ 
proximately 10-30 p,Jy. The radio through IR/optical to 
X-ray spectra of these three regions of the 3C 31 jet are 
very different, but each region mirrors in some manner 
features seen in oth er jets le. g. M87 (IShi et al. 


iForman et al. 11200711. Cen A ( Hardcastle et al. 
3C 66B ( Tanslev et al. I I^OOOL i Hardcastle et al. 




2006), 


2001 ), 


and 3C 15 (iMartel et a/. 11199^ 1. The change in the jet 
from X-ray dominated to IR dominated further down 
the j et is also seen in jets of greater power (e.g. 3C 
273; lUchivama et al. l[2006f) . In the region closest to the 
nucleus, Rl, the X-ray flux is the highest, the break fre¬ 
quency is the lowest, and the spectrum is the flattest after 
the break. The relatively flat spectral index extending to 
the X-ray band is similar to tha t seen in the knots of th e 
inner kpc of the jet in Cen A (|Hardcastle et aL]l2003D . 
Further extending the analogy to the jet of Cen A, it is 
likely that this region of the 3C 31 jet contains significant 
unresolved substructure in the X-ray and probably the 
optical/IR as well. 

The IR emission in regions R2 and R3 is significantly 
enhanced relative to the X-ray in comparison with re¬ 
gion Rl. The break frequencies are also higher than in 
Rl (see Tabled, as are the post-break spectral indices. 
The simplest explanation for the observed radio, IR, op¬ 
tical, and X-ray fluxes in the three regions is that there 
is a single impulsive event in region Rl, and the spectral 
steepening is simply the result of progressive radiative 
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losses. The radiative lifetime of the relativistic particles 
emitting in the Spitzer band is ^3x10^ yrs, assuming 
they are radiating in the equipartition field. If these par¬ 
ticles are advected down the jet with velocity 0.5c, they 
would travel roughly 5 kpc. The IR emission from the jet 
is detected to at least 10 kpc from the nucleus, so based 
simply on this lifetime argument, it is plausible that the 
observed multi-frequency features of the jet are the result 
of synchrotron aging. However, the X-ray and IR mor¬ 
phology of the jet argues strongly against this scenario. 
In particular, the jet is clearly brightest in the IR near 
the boundary between regions R2 and R3, and fainter 
closer to the nucleus within region R2. This suggests 
that (re-)acceleration is occurring near the boundary of 
R2/R3. The IR flux in region RI is not well constrained 
(we only infer it from the radio/optical/X-ray spectrum 
due to contamination from the nucleus), but the pro¬ 
gressive brightening of the IR emission with increasing 
distance from the nucleus demonstrates that the jet is 
being re-energized in regions R2/R3. Additionally, there 
would be no X-ray synchrotron emission in region R2 in 
this model. Given the radiative lifetimes, advection of 
ultra-relativisitic particles may play some role in the ap¬ 
pearance of the jet, but it clearly is not the whole story. 

There must be a second site of particle acceleration in 
regions R2/R3. The dramatic change in spectral shape 
between Rl and R2 suggests that the jet is much less ef¬ 
ficient at accelerating particles to 7 ^I 0 ^“® perhaps be¬ 
cause it has significantly decelerated. Additionally, the 
morphological changes in the radio also suggest a sig¬ 
nificant change in the flow dynamics from region Rl to 
R2/R3. The calculations of L02 show a sharp drop in 
velocity between regions Rl and R2, and a more gradual 
deceleration of the jet between R2 and R3. The obvious 
question to ask is are the changes in particle acceleration 
efficiency related to the morphological/dynamic changes 
in the jet and if so, how are they related. More gener¬ 
ally, the broad-band spectral features seen in 3C 31 are 
similar to that seen in a wide variety of radio galaxies. 
The spectra of other extragalactic jets generally show a 
low frequency spectral index between 0.5 and 0.7 and 
are typically fit with a broken power law. The jets of 
3C 66 B, 3C 15, and the outer diffuse regions of the Cen 
A jet all have breaks in their spectra around 10^^ Hz 
and have high frequency spectral indices between 1.13 
(Cen A, 3C 15) and 1.35 (3C 66 B), similar to region R2. 
We speculate that the efficiency of ultra-relativistic par¬ 
ticle acceleration is a strong function of the jet velocity 
and/or the entrained mass, and that this is a general 
feature of particle acceleration of all extragalactic jets 
on kpc scales. 

Mass loading may have an important effect on the dy¬ 
namics of particle acceleration. In both Cen A and 3C 31, 
the jets are believed to be moderately relativistic over the 
region where X-ray and IR synchrotron emission is de¬ 
tected. In both cases, the mass swept up by the jet, either 
from embedded AGB stars or from the ambient ICM, is 
relatively small. But in both cases, the IR synchrotron 
emission extends considerably further from the nucleus 
than the X-ray emission. There are many ways that a 
small amount of mass loading can fundamentally change 
the dynamics of the flow: it could change the internal 
sound speed, the equation of state, the Reynolds num¬ 
ber, etc. A full discussion of the effects of mass loading 


the jet on the particle acceleration is beyond the scope 
of this paper, but given the similarities of 3C 31 and Cen 
A, the dynamics of these systems appear to be qualita¬ 
tively similar. More generally, deep radio observations 
of FR I jets typically show knotted embedded structures 
(|Worrall et al. l[2007ll similar to Cen A and 3C 31, sug¬ 
gesting that dissipative processes are a common feature 
of FR I jets and play an important role in particle accel¬ 
eration. 

The dramatic transition from the X-ray dominated to 
IR dominated jet is most clear in Centaurus A where the 
X-ray jet narrows and disappears roughly 4 kpc from the 
nucleus, while the IR emissio n brightens as the jet op ens 
into the northeast radio lobe (iHardcastle et al. fcoodl . It 
is in this region where the jet dynamics appears to transi- 
tion from momentu m-dominated to pressure-dominated. 
IKraft et al. I (120081 1 showed that in the case of Cen A, the 
jet appears to cross a density/temperature discontinuity 
in the hot ISM. This discontinuity in density and temper¬ 
ature, but not pressure, is the result of non-hydrostatic 
motions induced in the gas by a recent merger. There 
are several galaxies in the larger 3C 31 group, and one 
of these satellites may have disturbed the gas core. The 
transition from momentum to pressure domination in the 
jet may have an important effect on the particle acceler¬ 
ation efficiency and hence the visible appearance of the 
jet at IR/optical/X-ray wavelengths. This suggests that 
subtle changes in the state of the external atmosphere 
may have an important influence on the efficiency of elec¬ 
tron acceleration internal to the jet. The existing Chan¬ 
dra data are not sufficient to determine whether there is 
a surface brightness discontinuity present in the gas at 
the interface between regions Rl and R2. Alternatively, 
it has been argued that the X-ray synchrotron shocks in 
the Cen A jet are the result of interactions between the 
jet and AGB stars in the host galaxy (Nulsen et al., in 
preparation), and that mass loading of the jet from mate¬ 
rial swept up via Kelvin-Helmholtz instabilities plays lit¬ 
tle role in the deceleration of the jet. There are some im¬ 
porta nt differences between the Cen A jet (jKraft et al. I 
l2002fl and that of 3C 31; the jet is somewhat smaller in 
Cen A and its synchrotron output is about an order of 
magnitude less powerful than in 3C 31, so one must be 
careful about pushing this analogy too far, but there are 
clearly some qualitative similarities. 

We obtained the external gas pressure measured 
as a function of di stance from the nucleus from 
iLaing fc Bridld (|2002a[l , and compared it to the equipar¬ 
tition pressure for each region. In region Rl, the dif¬ 
ference between the external pressure and the equipar¬ 
tition pressure is nearly an order of magnitude, but in 
region R3, it is less than a factor of 2. Note that the 
jet pressures cannot actually be this low, the jet would 
be crushed in an (IS M) sound crossi n g time if this were 
true. As argued bv iLaing &: Bridld (l2002all . the trend 
that Peq/Pgas approaches unity at larger distances from 
the nucleus suggests that mass entrainment and heating 
is not providing the additional pressure support. Rela¬ 
tivistic, radiating particles and magnetic fields are likely 
the dominant source of pressure support along the length 
of the jet. One would naively expect that more mass has 
been entrained at larger radii and that if such entrain¬ 
ment played an important role in the pressure of the jet, 
it would become progressively more important at larger 
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radii. That is, mass loading cannot play a significant 
role in providing pressure support, although as we have 
argued above it likely plays a critical role in particle ac¬ 
celeration. 

Finally we note that the clear differences in radio flux 
density between the jet and counterjet can be acco unted 
for by relativistic beaming (|Laing fc Bridlell2002all . and 
the differences in the IR can be attributed to the same 
process. We calculated the expected flux ratio in the jet 
and counterjet due to beaming: 


/ I -/3 

\1 -I- /3 cos{9)) 


( 1 ) 


Assuming an angle to the line of sight of 6* = 52° and 
velocities from the models of L02 and using the appro¬ 
priate spectral indices, we find the expected ratio in the 
region where we measured the counterjet flux to be ^3 in 
the radio bands and ~6 in the infrared bands. Our mea¬ 
sured ratios are close to these in the radio bands and the 
marginal detections/upper limits on the counterjet’s in¬ 
frared flux yield lower limits on the flux ratio that agrees 
with the expected value. Since the X-ray emission is 
strongest in Rl, we searched for X-ray counterjet emis¬ 
sion in a corresponding region on the opposite side of the 
nucleus. We do not currently detect the counterjet. Us¬ 
ing the same velocity profile (/3 between 0.72-0.76) and 
angle to the line of sight, we expect a ratio between 50 
and 65 estimated over Rl. The expected (0.5-2.0 keV) 
flux of 0.1 njy is lower than the upper limit of 0.6 nJy 
currently set by the background in the existing Chandra 
exposure. 


5. CONCLUSIONS 

We have presented a clear detection of IR emission 
from the northern jet and a marginal detection of the 
counter jet of 3C 31. We show that there is significant 
variation in the radio through X-ray spectra along the 
jet, and have outlined some implications for the jet dy¬ 


namics and particle acceleration. In particular, the ratio 
of IR to X-ray flux increases dramatically as a function 
of distance from the nucleus. In a future paper, we will 
present results from Spitzer observations of a sample of 
nearby radio jets (Bliss et al, in preparation). We are 
slowly increasing our sample size of IR, optical, and X- 
ray jets, and it is now clear that broadband synchrotron 
emission is a ubiquitous feature of FR1 jets. The spectral 
properties of the 3C 31 jet appear to span those seen in 
other extragalactic jets. It is likely that a range of parti¬ 
cle acceleration regimes are present in jets generally, and 
at least two different regimes are present in 3C 31. 

Additional future observations may shed further light 
on some of the issues addressed here. In particular, a 
deep Chandra observation would be able to resolve the 
detailed structure of the gas external to the jet. Such an 
observation may reveal contact discontinuities between 
two relatively stationary fluids that may play an impor¬ 
tant role in the jet dynamics. A deeper X-ray observation 
of the jet may also reveal regions of lower surface bright¬ 
ness further from the nucleus than observed here, as well 
as more clearly resolve the structures of the inner jet. 
Such an observation also would put strong constraints 
on the presence of a counter jet and constrain the jet 
velocity. Higher angular resolution IR observations, par¬ 
ticularly at lower frequencies than presented here, may 
be able to resolve knots and other substructures in the 
jet. 


This work was based on archival data obtained from 
the Spitzer Science Archive, the Chandra Data Archive, 
and the Multimission Archive at STScI. This work was 
supported in part by the Smithsonian Institution, the 
Chandra X-ray Center, NASA contract NNX07AQ18G, 
and the Harvard College Observatory. 

Facilities: Spitzer, CXO, VLA, HST:WFPC2, 2MASS, 
GALEX 
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